Population Inversion Induced by Resonant States in Semiconductors 
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We present a theoretical prediction of a new mechanism for carrier population inversion in semi- 
conductors under an applied electric field. The mechanism is originated from a coherent capture- 
emission type inelastic scattering of resonant states. We support our theory with concrete calcula- 
tions for shallow acceptor resonant states in strained p-Ge where a lasing in THz frequency region 
has been recently observed. 
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Besides elastic scattering, resonant states also cause a 
specifically strong inelastic scattering because electrons 
can have a finite lifetime in the resonant state. Con- 
sequently, the electron energy is not conserved, and it 
can be re-emitted with the energy different from the ini- 
tial value. We will prove in this Letter that due to such 
an effect a novel population inversion of carrier distribu- 
tion can be realized in non-equilibrium semiconductors 
in the so-called streamer regime. The streamer regime 
emerges in a relatively pure semiconductors at low tem- 
peratures where both impurity scattering and acoustic 
phonon scattering are weak. Under an applied electric 
field a carrier drifts in the momentum space almost scat- 
ter free until its energy reaches the optical phonon energy 
Hluq. Then the carrier emits an optical phonon and re- 
turn to the low-energy region jy . If there exist a resonant 
state with energy Eq<Tiujq, it acts as a carrier trap with 
a finite lifetime which is field- independent. In the energy 
space one then finds carriers accumulation in the vicinity 
of the resonant energy Eq. 

We will use a general model to analyze this novel mech- 
anism for population inversion. While our theory applies 
to a variety of resonant states, in this Letter we will also 
demonstrate in quantitative details that our theory is 
realistic. For this purpose we will investigate the popula- 
tion inversion in uniaxially strained p-Ge where shallow 
acceptors induce resonant states f|,[|. It is important to 
notice that based on such a p-Ge system, a laser oper- 
ating in terahertz frequency region has been fabricated 
recently [Q. By studying the emission lines positions as 
functions of applied stress, it was confirmed m that 
the radiation emission is due to optical transitions be- 
tween the resonant states and the localized acceptor (Ga) 
states. However, the mechanism of lasing has remained 
a puzzle. Our theory will explain the physical origin of 
the population inversion which leads to the lasing. 

We consider a model system of charge carriers inter- 
acting with optical phonons of energy Hujq under a strong 
external electric field £ along the z axis. In this case an 
electron acquire energy from the electric field and drift in 



the k-space until its energy exceeds Hloq. Then the elec- 
tron emits an optical phonon and returns to the region 
fcwO. Thus, our model includes a drain D at E\^=%ujo and 
a source S at small k. In the absence of resonant scat- 
tering by impurities, for carrier kinetic energy E-^-^-huQ, 
the carrier distribution function /k can be obtained from 
the kinetic equation Q) 

dfk , e£ <9/ k 



dt 
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= S-D. 
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For our problem we can well approximate the drain 
D by a black-wall boundary condition that /k=0 for 
Ek>hujQ. The source intensity S is determined by the 
k-space particle flow with energy Ek=hu)Q. In the fol- 
lowing we use the expression for S as 



S = S (t)Q(e 
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with m z being the effective mass along z-axis. The fre- 
quency va is related to the rate va y {E^/Tlloq) — 1 of op- 
tical phonon emission by the carriers with E^hujQ. The 
source amplitude So(t) is determined from the condition 
of particle flow conservation in the k-space, namely the 
flow out of the surface E^=Tluiq at time t returns evenly 
back into the region _E"k<eoj 



e J(e-dS)f k (t) 
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The integration is performed over the surface defined by 
the equation Ek=hujQ. 

The stationary solution of Eq. (jl]) (the so-called 
streamer) is that /k is almost constant if k lies in a cylin- 
der, and /k=0 otherwise. This cylinder in k-space is 
determined by 



< k z < 
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where m± is the transverse component of the effective 
mass tensor. 

When impurities induce resonant states with complex 
energy E +iY/2, charge carriers may be trapped at en- 
ergy E <Twjo for a time interval ~ U/T. As a result, a 
maximum of the non-equilibrium distribution function is 
formed around the energy Eq. This is just the population 
inversion. To take into account the resonant scattering 
by impurities of concentration Ni , at the right hand side 
of Eq. (Q), we should add the impurity collision integral 
I, 



I = NiV^2 [/k-WW ~ /kW k - k ] 

k' 
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where V is a normalization volume, and f r is the resonant 
state population which satisfies the kinetic equation 

dfr = \ - 
dt ~ ^ 

k' 

The first term in Eq. (||) represents the elastic scattering, 
while the second one describes the coherent capture and 
re-emission by the resonate state. The sum W r k> 
is the total escape rate from the resonant state, and is 
equal to Y/fr. 

The transition probabilities Wkk' and Wk r can be ex- 
pressed through the respective scattering amplitude ikk' 
and transition amplitude tkr as 
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tkk' an< i tk r can be calculated using the Dirac ap- 
proach H to the scattering problem at a resonant state. 
According to this approach, via hybridization with ex- 
tended states, the wave function <p(r) of the bare local- 
ized impurity state develops into the resonant state. The 
general form of the scattering state is 



*k(r) = -j=e 



kr 



^kk' 
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Finally, we will add the normalization condition 

£ /k + NiVfr = nV (10) 

k 

where n is the total electron concentration, and solve all 
these equations selfconsistently. 

The above theory is general for different types of reso- 
nant states. To demonstrate explicitly the population in- 
version predicted by our theory, we will investigate quan- 
titatively the resonant states induced by shallow accep- 
tors in uniaxially strained p-Ge. 



In cubic semiconductors with symmetry group Oh the 
shallow acceptor wave functions 0^ M '(r) are 4-fold de- 
generate with the total angular momentum projections 
M=±l/2 and ±3/2. The ground state has Y% symmetry. 
Under a uniaxial strain the valence band top splits into 
two doubly degenerate energy levels. The correspond- 
ing two sets of wave functions transform according to Y^ 
and Ty" representations if the stress is along the [001] di- 
rection, and according to Y^ and Ye~, 6 representations 
if the stress is along the [111] direction. The ground 
state acceptor wave functions split in the same way and 
can be classified by the total angular moment projections 
(M=±l/2 or M=±3/2) along the stress direction. In our 
calculations, we will use the spherical approximation for 
the Luttinger Hamiltonian (LH) |(| . In the limit of large 
strain such that the splitting Ed at the top of subband 
exceeds the Coulomb energy, the LH can be treated with 
a quasi-diagonal approximation, and is represented by 
two 2x2 blocks. The states in each subband can now be 
classified by the projections (to = ±1/2 or ±3/2) of the 
hole spin on the stress axis fm . 

In this way, in each of doubly-degenerate subbands, we 
obtain both extended states and the localized Coulomb 
states below the subband. The corresponding energy lev- 
els and wave functions are calculated following the vari- 
ational procedure in Ref. |^||. For a large strain, the 
energy levels are shown in Fig. [I]. 
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FIG. 1. The acceptor levels diagram for uniaxially strained 
Ge in the large strain limit along [001] . The optical transitions 
observed in the lasing U are indicated by arrows. 



The off-diagonal terms of LH mix the states belong- 
ing to different subbands. Treating these terms as per- 
turbation in the Dirac approach Q] the hybridization 
of local and extended states yields the resonant state 
with energy E and level width Y. Along the Dirac 
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approach one also obtains the scattering amplitudes 
t^g? (m, m'=±l/2) and the transition amplitudes i^ l r m 
(to=±3/2 and m'=±l/2). The explicit expressions of 
these amplitudes, as well as the details of deriving them 
will be published elsewhere. 

We are ready to use these scattering and transition 
amplitudes to solve the set of kinetic equations (0)-(||)- 
However, for p-Ge the kinetic equations can be simplified. 
Let T£ be the transient time during which the stationary 
non-equilibrium distribution is established. Tg depends 
on both the electric field £ and the impurity concentra- 
tion iVj. As will be shown below, for p-Ge, there exists 
rather large region of electric fields and impurity con- 
centrations where rg is much larger than the lifetime of 
the resonant state Ti/T. In this case we can set the left 
hand side of Eq. (^) to zero, and so the occupation of 
the quasi-local states, f r , follows adiabatically the dis- 
tribution function /k of the extended states. Since the 
localized and the extended states are doubly degenerate, 
we have / r f3/2 =/ r T 3/2 =/ r and f£ 1/2 =/ k 1/2 =/k- If we 
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we obtain 
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Substituting Eq. @ into Eqs. (0) and (§), we arrive at 
the kinetic equation for /k 
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where 1 1. 
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The boundary 



condition for the above kinetic equation is / k =0 at 
E^=fiujQ. The source So(t) and the energy eo are given 
by Eqs. (|]) and (j|), respectively. 

We will solve Eq. (|l^) numerically as a non-stationary 
equation. We start from some initial distribution and 
then follow the evolution of the distribution function until 
it reaches a stationary one. In this way, the final station- 
ary distribution and the transient time rg are obtained 
for various values of applied electric field and impurity 
concentration. In our numerical calculation, the values 
of material parameters for p-Ge are Luttinger parame- 
ters 7i=13.38, 72=4.24, 73=5. 69, hu — 36 meV, and 
the characteristic frequency zm=5x10 12 s _1 . 
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FIG. 2. Normalized distribution function (a/2ir) 3 n~ 1 
as a function of squared normalized momentum (k z /a) 2 (up- 
per panel) and (k±/a) 2 (lower panel) in a uniaxially strained 
p-Ge under an electric field £ = 100 V/cm along the z axis. 
Various curves are for impurity concentration Ni=l (curve 1), 
1.1 (curve 2), 1.2 (curve 3), 1.4 (curve 4), and 1.5 (curve 1), 
in units 10 15 cm" 3 . 



We will set the pressure at 5 kbar and the electric field 
at £=100 V/cm along [111]. In this case, the resonant 
level is at the energy Eo = 10 meV and has a width 
r = 2 meV. Here we have adopted the approximation of 
a single hole band, and we are aware of the fact that at 
a pressure of 5 kbar the inter-hole-band splitting £^=20 
meV is less than Hujo- One can show that including the 
second band will only change the tail of the distribution, 
and such change is not important for the problem under 
consideration. From Eq (0) we obtain the source width 



£o=4.25 meV. If we define a—y 2nioEo/h 2 ji as the unit 
for wave vector, where mo is the free electron mass, our 
calculated normalized distribution functions are shown in 
Fig. ^| as functions of normalized k z /a (upper panel) and 
of normalized k±/a (lower panel). The curves in Fig. 2 
are for impurity concentration Ni=l (curve 1), 1.1 (curve 
2), 1.2 (curve 3), 1.4 (curve 4), and 1.5 (curve 1), in units 
10 15 cm~ 3 . While the distribution function peaks are 
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centered at the resonant energy Eq, their corresponding 
peak positions in the upper panel are different from those 
in the lower panel. This is due the anisotropy of the effec- 
tive mass: m z =0.04mo and m±=0.13mo. For the results 
shown in Fig. |2|, the transient time is about T£«fO~ n s, 
which is much longer than the lifetime ?i/r«2xfCP 13 s. 
Consequently, our calculation based on the condition that 
Tre^h is selfconsistcnt. 
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FIG. 3. Normalized distribution function as a function of 
squared normalized momentum (k z /a) 2 for impurity concen- 
tration Ni — 5 x fO 15 cm -3 . The curves are: (a) no resonant 
scattering, (b) include only elastic scattering from resonant 
states, and (c) include full resonant scattering. 

In order to demonstrate precisely which scattering pro- 
cess is responsible to the population inversion, we have 
used the same values of material parameters and same 
electrical field strength to calculate the distribution func- 
tions for 3 cases, and the results are shown in Fig. ^. In 
the absence of resonant states, the distribution function 
follows curve (a). When elastic scattering due to reso- 
nant states is taken into account, the result change into 
the step-like curve (b). Finally, by adding the inelas- 
tic scattering process, a peak emerges in curve (c). It 
is then clear that the population inversion is induced by 
the capture- emission component. The population accu- 
mulated in resonant states is controlled by the distribu- 
tion function in the continuous spectrum, as indicated by 
Eq. (fi"f|). At the same time, the non-equilibrium popu- 
lation in the localized states can be only less than that 
in the low-energy continuous states. Consequently, un- 
der the conditions imposed on our calculation, the intra- 
center population is also inverted. It has been demon- 
strated in Ref. @|| that the lasing in uniaxially strained 
Ge:Ga is connected to the transitions shown by arrows 
in Fig. [j]. Among these transitions, besides the allowed 
intra-center optical transitions between the lowest res- 



onant s-type state and excited localized p-type states, 
there is also transition corresponding to the forbidden s- 
type to s-type. The observed forbidden transition is due 
to the accumulation of the continuous spectrum carriers 
at the energy around Eq. These states are almost plane 
waves and so turn the forbidden process into an allowed 
one. The sharp emission line in the observed optical spec- 
trum proves the existence of the peak in the distribution 
function /k. Thus, the resonant-states-induced popula- 
tion inversion predicted by our theory explains the origin 
of lasing in terahertz frequency range observed in these 
strained Ge:Ga samples. 

In conclusion, we have predicted that resonant states 
can produce a population inversion in the carrier dis- 
tribution function in strained semiconductors under an 
external electric field. Our theoretical prediction is con- 
firmed by concrete calculations for strained p-Ge, where 
resonant states give rise to the lasing observed in THz fre- 
quency region. We believe that the proposed mechanism 
for population inversion is rather general since resonant 
states can be created by various means. 
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